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variability  in  the  upper  layers  of 

the  central  western  Mediterranean,  MILEX-82: 

Final  report 


J.  Moen 


r' 

Abstract:  The  1982  Mixed-Layer  Experiment  (MlLEX-82fjis  described 
in  detail  and  the  data  are  presented  in  ways  designed  to  illustrate  the  spa¬ 
tial  and  temporal  vaii^ility  of  the  surface  and  deeper  layers  of  the  central 
western  Mediterranean^uiing  the  autumn  transition  period.  The  surface 
mixed-layer  is  seen  to  coot  and  to  deepen  due  to  ronvection  and  mechanical 
stirring  induced  by  passing  storms.  . 

-The  predominating  motions  observed;  and  the  ones  given  most  attention 
In  this  report,  are  the  near  inertial  internal  waves  propagating  horizontally 
and  vertically  from  the  mixed  layer  into  the  pycnocline.  Possible  genera¬ 
tion  mechanisms  for  the  vertical  oscillations  are  discussed,  it  is  argued  that 
these  may  originate  weU  to  the  north  of  the  study  are^in  the  Golfe  du  Lion,, 
where  a  high  amplitude,  inertially  oscillating,  Ekman  mass-flux  divergence 
can  be  expected  due  to  the  high  values  of  wind-stress  curl  in  that  region. 

Data  are,  in  some  cases,  faith'*  extensively  hvorked  up*  to  produce  a  vari¬ 
ety  of  contour  plots  for  temperature,  salinity,  density  and  dynamic  anomaly 
and  a  variety'oPcomposite  time-series  for  both  the  thermistor-chain  data 
and  current-meter  array  vectors.  Also  presented,  but  discussed  in  less  de¬ 
tail,  are  a  number  of  measured  and  derived  quantities  relevant  to  air-sea 
interaction  and  mixed  layer  deepening.  It  is  considered  that  these  will  be 
useful  for  comparison  with  other,  similar  experiments  and  with  output  from 
numerical  models.  *r  , . 
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1.  Introduction 


A  mixed-layer  experiment  (hereafter  referred  to  as  MILEX-82),  designed  to  obtain 
a  three-month  coverage  of  the  surface  layer  structure  (upper  100  m)  at  one  location, 
was  begun  in  late  September  1982  in  the  central  western  Mediterranean.  Due  to 
a  combination  of  weather,  mooring  technique  and  instrumentation  used  (to  be  dis¬ 
cussed  later)  the  experiment  had  to  be  terminated  earlier  than  planned,  shortening 
the  duration  to  just  over  one  month.  Fortunately,  this  was  still  a  useful  timespan 
for  a  mixed-layer  experiment  -  since  considerable  mixing  and  deepening  activity 
occurred  throughout  the  period. 

The  purpose  of  MILEX-82  was  threefold;  firstly,  to  investigate  the  dominant  pro¬ 
cesses  responsible  for  surface  layer  mixing  in  the  Mediterranean;  secondly,  to  obtmn  a 
good  data  set  against  which  to  test  and  adapt  one- dimensional  numerical  mixed-layer 
models;  and  thirdly,  to  complement  the  data  gathered  during  the  earlier  MILEX-80 
experiment  in  the  north  Tyrrhenian  Sea. 

It  became  evident  during  analysis  of  MILEX-80  data,  described  in  Moen  (1984), 
that  one  fairly  important  objective  could  not  be  adequately  fulfilled:  despite  the 
overall  quality  of  the  data  and  the  general  success  of  the  experiment,  the  complexity 
of  the  local  oceanography  made  the  MILEX-80  region  unsuitable  for  the  testing  of 
one-dimensional  mixed-layer  models.  Upwelling  and  horizontal  advection  tended  to 
dominate  vertical  mixing  in  determining  the  depth  of  the  mixed  layer.  Although 
gross  estimates  of  upwelling  were  made  from  the  data,  they  could  not  be  separated 
from  those  of  vertical  mixing  to  the  resolution  necessary  for  one-dimensional  model 
testing. 

Thus  it  was  hoped  that  MILEX-82  would  provide  data  in  which  horizonteil  advec¬ 
tion  and  upwelling  were  either  effectively  absent  or,  if  present,  completely  resolvable 
due  to  proper  choice  of  space  and  time  sampling.  Furthermore,  the  data  were  to 
be  obtained  under  environmental  conditions  similar  to  those  of  MILEX-80,  i.e.  dur¬ 
ing  the  same  season,  and  during  similar  air  and  sea  temperattires,  similar  wind 
speeds,  humidities,  cloud  covers,  radiant  heating  and  cooling  conditions,  and  similar 
surface-layer  stratifications.  In  this  way  a  mixed-layer  model  ‘tuned’  to  MILEX-82 
conditions  might  be  applied  directly  to  MILEX-80  data  for  the  purpose  of  separating 
the  purely  mixing  influences  on  the  mixed-layer  depth  from  the  purely  dynamic. 

It  is  important,  from  the  military  oceanography  viewpoint,  to  recognize  that  the 
dynamics,  in  the  form  of  upwelling  and  advection,  are  often  more  effective  than 
mixing  in  determining  the  depth  of  the  mixed  layer.  This  is  indeed  the  case  for 
much  of  the  Mediterranean  where  the  wind  stress  is  known  to  have  a  high  spatial 
variability  (as  an  example  see  Moen,  1934). 
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Another  limitation  of  the  MILEX-80  data,  from  the  point  of  view  of  mixed-layer 
modelling,  was  its  rather  short  duration  of  approximately  three  weeks.  The  longer 
time  series  of  MILEX-82  will  not  only  be  more  effective  for  model  testing  (not 
included  in  this  report )  but  will  also  better  resolve  the  influences  of  lower  frequency 
effects  on  the  surface  layer.  These  effects  had  periods  ranging  from  the  inertial 
(18.66  h)  to  the  passing  of  cyclonic  lows  (several  days). 

MILEX-82,  then,  was  performed  with  the  combined  intentions  of  extending  and 
complementing  the  data  base  of  MILEX-80.  It  was  conducted  using  mainly  existing 
instrumentation  and  with  the  benefit  of  data  preparation  and  analysis  programs 
previously  developed  for  MILEX-80. 

The  report  therefore  presents  analyses  of  the  MILEX-82  data  and  indicates  their 
potential  for  comparison  with  mixed-layer  models  in  a  Mediterranean  environment . 
Much  of  the  interpretation  of  the  data  concerns  the  properties  and  effects  of  inertial 
oscillations  in  the  mixed  layer  and  thermocline. 
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2.  Description  of  experiment 


2.1.  DEPLOYMENT  OF  INSTRUMENTS 

Three  thermistor-chain/current-meter  moorings,  marked  Al,  B1  and  Cl  in  Fig.  1, 
were  deployed  in  a  triangular  configuration  of  side  10  km  with  centre  at  40°00'N, 
006°00'E  in  the  central  western  Mediterranean.  The  separation  was  chosen  to  corre¬ 
spond  to  the  local  Rossby  radius  of  deformation  and  hence  the  expected  length  scale 
of  variability  at  this  location.  Three  additional  current-meter  moorings,  marked  A2, 
B2  and  C2  in  Fig.  1,  were  deployed  within  close  proximity  (<  1  km)  of  Al,  B1 
and  Cl  respectively.  All  moorings,  shown  schematically  in  Fig.  2,  were  of  the  taut- 
nylon/surface-float  type  with  the  instruments  arrtinged  to  span  the  surface  mixed- 
layer  and  thermocline.  Aanderaa  meteorological  chains  were  attached  at  3.5  m  to 
the  surface  floats  of  moorings  B1  and  Cl.  The  position  of  each  mooring  is  given  by 

Al;  39'56.00'N,  006“31.70'E 
A2:  39’56.60'N,  006“31.90'E 
Bl:  40*00.01'N,  006"26.15'E 
B2;  40“00.30'N,  006‘'25.20'E 
Cl:  40*02.10'N,  006'32.70'E 
C2;  40'02.20'N,  006'32.90'E 

The  choice  of  spatial  distribution  for  the  moorings  and  instruments  was  influenced 
by  two  factors.  Firstly,  should  they  become  important,  first  spatial  derivatives  at 
several  different  levels  could  be  derived  and  hence  crude  estimates  made  there  of 
horizontal  advection  and  divergence.  Secondly,  should  they  not  be  important  (i.e.  in 
the  absence  of  advection,  divergence  and  strong  horizontal  gradients  -  which  turned 
out  to  be  the  case),  the  moorings  would  be  sufficiently  close  that  horizontal  separa¬ 
tion  could  be  ignored  and  greater  vertical  resolution  obtained.  In  particular,  such 
parameters  as  vertical  shear,  bulk  gradient  -  Richardson  number,  and  frequency, 
phase  and  energy  of  inertial  oscillations  could  better  be  studied  as  a  function  of 
depth.  Another  advantage  of  the  mooring  configuration  was  in  its  duplication  of 
measurements,  particularly  with  regard  to  the  thermistor  and  meteo-chains.  This 
turned  out  to  be  well  justified  since  three  of  the  moorings  came  adrift  in  severe 
weather  and  the  experiment  had  to  be  terminated  much  earlier  than  planned  (see 
Sect.  2.2.). 
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2.2.  DURATION  OF  EXPERIMENT 

The  moorings  were  deployed  during  the  last  week  of  September  1982  and  recovered 
during  the  first  week  of  November,  with  the  longest  records  being  of  approximately 
40  days  duration  Figure  3  presents  a  timetable  for  each  of  the  instruments,  includ¬ 
ing  times  of  XBT  and  CTD  casts.  Also  shown  in  Fig.  3  are  depths  and  sampling 
intervals  for  each  of  the  moored  instruments.  In  the  case  of  the  thermistor  chains 
the  depths  refer  to  the  uppermost  and  lowermost  thermistors.  Chains  of  22  thermis¬ 
tors  were  formed  by  joining  two  11-element  Aanderaa  thermistor  chains  end-to-end. 
The  spacing  allowed  highest  vertical  resolution  in  the  vicinity  of  what  was  guessed 
would  be  the  niixed-layer/thermocline  interface  region  as  the  mixed  layer  deepened. 
The  depths  of  the  thermistors  are  thus  given  relative  to  the  uppermost  progressing 
downwards,  and  relative  to  the  lowermost  progressing  upwards,  by  the  sequence  9.7, 
12.7,  15.7,  18.7,  21.7,  24.7,  27.7,  32.7,  37.7  and  42.7  m. 

Attached  in-line  to  five  of  the  moorings  (see  Fig.  2)  were  NBA  current  meters.  This 
turned  out  to  be  an  unfortunate  oversight  in  the  mooring  design  and  the  main  reason 
for  terminating  the  experiment  early.  It  was  suspected  from  measurements  made 
during  MILEX-80  that  the  NBA  instrument  might  quite  seriously  over-estimate  cur¬ 
rent  speed  when  attached  to  a  taut-nylon  surface-mooring,  but  it  was  not  suspected 
that  the  vertical  oscillations  of  the  mooring  line  could  be  so  energetic  as  tc  snap 
the  tension/tie  bar  of  the  current  meter.  It  appears  that  (due  to  large  direction 
vanes  on  its  tail  end)  asymmetry  of  the  drag  forces  on  the  current  meter  forced 
it  to  oscillate  vertically,  causing  an  upward  or  downward  tilting  of  the  instrument 
that  pointed  always  into  its  direction  of  motion.  Since  this  periodic  tilting  could  be 
arrested  only  by  the  tension/tie  bar  passing  through  the  centre  of  buoyancy  of  the 
current  instrument,  a  critical  stress  occurred  there  -  causing  the  bar  to  fracture  and 
snap.  This  stress  was  related  to  the  energy  and  d"ration  of  the  vertical  oscillations, 
which  in  turn  were  related  to  the  sea  state  and  possibly  were  augmented  through 
resonance  of  the  stretched  nylon  rope  with  one  of  the  lower  frequency  surface  waves. 

MILEX-82  took  place  in  three  phases:  a  mooring  deployment  phase,  conducted  by 
the  SACLANTCEN  research  vessel  Maria  Paolina  C  during  which  XBT  and  CTD 
measurements  were  made  on  a  grid  as  indicated  in  Fig.  1;  a  CTD/XBT  survey  phase 
(Fig.  4a)  conducted  by  the  Italian  Navy  hydrographic  research  vessel  Ammiraglio 
Magnaghi,  during  which  surveillance  of  the  moorings  was  carried  out;  and  an  emer¬ 
gency  recovery  phase,  again  using  the  Maria  Paolina  G.  It  was  indeed  fortunate 
and  invaluable  to  the  experiment  that  the  Italian  Navy  were  in  the  vicinity  during 
the  period  when  moorings  broke  loose.  This  incident  illustrates  the  importance  of 
regular  surveillance  of  surface  moorings  in  open  ocean  experiments. 

Recovery  of  the  drifting  moorings  was  greatly  expedited  through  the  assistance  of 
the  Italian  Navy  aircraft  Spada  I  (30th  Air  Squadron)  in  response  to  a  request 
from  SACLANTCEN  to  MARISTAT  Rome.  The  ship’s  track  for  recovery  of  the 
moorings  is  indicated  in  Fig.  4b. 
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3.3.  DATA  RECORDING  AND  PREPARATION 


Sampling  intervals  for  each  of  the  instruments  arc  given  in  Fig,  3.  These  were  15  min 
for  the  vector  averaging  current  meters  (VACMs)  and  the  thermistor  chains,  6  min 
for  the  two  surviving  NBA  current  meters  (NBA566  and  567)  and  4  min  for  the 
two  Niel-Brown  acoustic  current  meters  and  for  the  vector  measuring  current  meter 
(VMCM). 

For  the  current  meters,  samples  represent  averages  over  the  sampling  interval,  thereby 
effectively  filtering  out  frequencies  greater  than  0.5t,  where  t  is  the  sampling  inter¬ 
val.  This  is  not  so  for  the  ther*nistor  chains:  as  supplied  by  the  manufacturer, 
they  have  a  time-constant  of  1  min,  and  so  sampling  at  less  than  1  cycle/min  could 
cause  aliasing.  This  problem  was  discussed  in  more  detail  in  Moen  (1984),  where 
it  was  demonstrated  that  ali'  ,ing  should  not  be  serious  providing  sampling  is  at 
If  cycle/min  or  higher.  In  order  to  increase  the  thermistor  time-constant,  thermis¬ 
tors  were  wrapped  using  a  4.5  mm  layer  of  PVC  tubing  and  tape.  The  frequency 
response  was  tested  by  applying  an  impulse  of  finite  width,  which  entailed  (540  s) 
plunging  the  thermistor  between  two  baths  at  different  temperature  and  measuring 
the  impulse  response. 


The  response  time  series  was  digitised  <md  Fourier  transformed  and  the  frequency 
response  function  H  was  calculated  by  the  formula 


«(/)  = 


I  g.y(/)  I 
I  GAf)  I  ’ 


where  G,^{f)  is  the  cross  spectrum  of  the  impulse  x  with  the  response  y  and  G,(f) 
is  the  autospectrum  of  the  impulse.  Figure  5  is  a  plot  of  H  which  shows  that  at 
4  cycles/h,  the  sampling  frequency  used,  the  thermistor  response  is  less  than  20%. 
Aliasing  for  the  thermistor  chain  data  was  therefore  not  a  problem. 


Data  from  the  moored  instruments  were  recorded  on  magnetic  cassettes,  transcribed 
onto  8-track  HP  tapes,  transferred  to  UNIVAC  disc  storage  and  then  prepared  for 
analysis  by  replacement  of  ‘bad’  data  values  and  trimcation.  Onboard  acquisition 
of  XBT  and  CTD  data  was  done  by  the  SACLANTCEN  lOS  (Integrated  Oceano¬ 
graphic  System)  on  the  Hewlett-Packard  21MX  computer.  Data  were  stored  on 
8-tratk  HP  tapes  and  later  transferred  to  UNIVAC  disc  storage  where  they  were 
‘cleaned’  by  removal  of  spikes  and  ‘bad’  profiles. 
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3.  Spatial  description 


3.1.  UPPER  LAYER  (ABOVE  150  m) 

The  MILEX-82  area,  by  contrast  with  that  of  MILEX-80,  was  found  to  contain 
relatively  little  horizontal  variation  in  respect  of  thermodynamic  variables. 


m  3.1.1.  Temperature  sections 

South-north  and  west-east  temperature  sections  for  the  upper  150  m  are  presented 
in  Fig.  6.  They  show  that  there  is  no  large-scale  horizontal  structure  in  the  south- 
north  direction  but  in  the  west-east  direction  the  thermocline  rises  appreciably  as 
it  approaches  the  shelf  break  to  the  west  of  Sardinia.  The  23  “C  isotherm,  for 
example,  which  represents  the  mixed-layer/thermocline  transition  depth  near  the 
mooring  triangle,  rises  from  about  35  m,  changing  gradually  at  first,  and  then  more 
rapidly  until  it  reaches  about  12  m  near  the  shelf  break.  It  would  appear  from  the 
west-east  section,  therefore,  that  the  depth  of  the  potentially-mixed  layer  increases 
fairly  rapidly  from  0  to  30  m  going  westward  from  the  Sardinian  shelf.  Just  how 
rapidly  is  not  possible  to  say  due  to  the  poor  spatial  resolution  afforded  by  the 
two  most  eastward  CTD  stations  10  and  1  on  Fig.  6b.  It  might,  for  example,  be 
attributed  to  an  internal  near-surface  front  crossing  somewhere  between  stations  10 
and  1. 

Although  isotherms  are  bending  upward  towards  the  shelf,  indicative  of  upweUing, 
the  surface  water  there  is  actually  warmer  than  that  to  the  west.  Such  a  situation 
could  not  prevail  long,  however,  in  the  presence  of  strong  winds  -  since  in  this 
case  the  surface  water  would  cool  quickly  due  to  vertical  mixing.  Figure  6b  thus 
provides  an  intersting  example  of  when  satellite  infrared  imagery  would  not  reveal 
the  underlying  structure  or  dynamics.  UpweUing  would  not  be  suspected  in  a  zone 
where  the  surface  temperature  was  seen  to  increase!  Fortunately,  this  anomalously 
cold  sub-surface  water  is  far  from  the  mooting  sites,  and  so  any  associated  dynamic 
effect  or  movement  of  the  structure  ought  not  to  interfere  with  the  one-dimensional 
mixing  assumed  to  be  taking  place  there.  A  better  indication  of  the  horizontal 
extent  of  the  cold  water  is  given  by  the  20  m  horizontal  section  of  Fig.  8a  (discussed 
below). 


■  3.1.2.  Salinity  sections 

A  slightly  different  picture  emerges  from  the  salinity  transects  shown  in  Fig.  7.  The 
south-north  section  (Fig.  7a)  indicates  low  salinity  gradients  crossing  the  mooring 
region,  involving  CTD  stations  9,  6,  5  and  7,  but  a  rather  strong  gradient  develops 
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throughout  the  whole  150  m  depth  south  of  station  26.  As  will  be  Skrgued  later, 
this  most  likely  corresponds  to  the  northward  branching  loop  of  the  North  African 
current  as  it  flows  eastward  from  the  Alboran  Sea  through  the  Sardinia  Strait. 

The  west-east  section  (Fig.  7b)  shows  a  surface  gradient  in  salinity  passing  right 
through  the  MILEX-82  zone,  strongest  near  the  surface  (~  0.5  ppt/100  km)  and 
becoming  very  weak  by  the  50  m  level.  Such  salinity  changes  do  not,  however, 
trsuislate  to  very  strong  density  gradients,  and  since  the  variations  are  also  very 
much  a  surface  phenomena  there  are  no  important  dynamics  associated  with  them. 
One  possible  explanation  is  the  following.  Water  slowly  upwelling  off  the  northwest 
Sardinia  coast  throughout  the  summer  would  tend  to  cause  a  shallower  mixed  layer 
there.  Heating  and  evaporation  would,  in  the  absence  of  strong  wind  mixing,  tend 
to  concentrate  heat  and  salt  in  the  shallow  layer  more  than  in  the  deeper  mixed 
layer  to  the  west.  Warmer  saltier  waier  could  then  spread  slowly  south-westward 
within  the  mixed  layer  due  to  Ekman  transport. 


■  3.1.3  Horizontal  sections 

Comparing  horizontal  sections  of  temperature  (Fig.  8a)  with  StJinity  (Fig.  9a)  at 
the  20  m  level  it  is  clear  that  the  salinity  maximum  occurs  to  the  west  of  the 
temperature  minimum.  The  closed  contours  indicated  in  salinity,  therefore,  do  not 
imply  the  existence  of  a  gyre  or  eddy  in  that  region.  They  might  imply,  instead, 
that  the  higher  salinity,  heavier  surface  water  present  above  the  doming  isotherms 
is  simply  displaced  westward  and  downward,  where  in  horizontal  section  it  shows  up 
as  a  region  of  anomalously  high  salinity. 

Comparing  both  the  temperature  (Fig.  8a)  eind  salinity  (Fig.  9a)  maps  at  the  20  m 
level  with  potential  density,  contoured  in  the  map  of  Fig.  10a,  it  is  clear  that  at  this 
level  temperature  has  the  dominant  effect  on  density.  Any  important  dynamics  will 
be  controlled  far  more  by  temperatme  than  by  salinity  -  a  situation  not  abnormal 
for  the  Mediterranean  surface  layer.  At  the  200  m  level,  however,  potential  density 
(Fig.  10b)  depends  equsdly,  if  not  more  strongly,  on  sabnity.  Verticed  sections  of 
density  have  not  been  included  since  they  appear  virtually  identical  in  form  to  those 
of  temperature. 


■  3.1.4.  Geostropic  currents 

To  examine  at  the  dynamics  of  the  MILEX-82  region  we  consider  maps  of  dynamic 
anomaly  csdculated  relative  to  motion  at  1000  m.  Two  such  maps  are  presented: 
Fig.  11a  for  20  m  (effectively  surface)  level  and  Fig.  11b  for  the  200  m  level.  Included 
also,  for  comparison,  is  a  wintertime  surface  dynamic-topography  map  (Fig.  11c) 
reproduced  from  Ovchinnikov  (1966). 

It  should  be  noted  here  that  in  common  with  ail  hori-.ontal  sections  describing  the 
MILEX-82  area,  the  contours  suffer  from  an  inadequate  grid  of  CTD  stations  - 
particularly  towards  the  extremities  of  the  grid.  Consequently,  any  interpretation  of 
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the  dynamics  to  the  east  near  the  colder  water  is  difficult  and  somewhat  subjective. 
There  is  a  tendeiicy  for  the  contouring  package  to  close  contours  going  away  from 
the  grid  region. 

Although  the  MILEX-82  maps  represent  only  one  survey  taken  in  the  autumn  period, 
whereas  Ovchinnikov’s  map  represents  wintertime  during  the  1950s  and  was  based 
on  a  coarser  survey  grid  (0.5°  squares),  the  two  pictures  show  some  agreement. 
In  both  maps,  for  example,  part  of  the  eastward  flowing  North  African  current 
loops  to  the  north  entering,  the  MILEX-82  zone  at  the  southwest  corner,  turning 
southeastward  and  then  leaving  via  the  southeast  corner.  The  MILEX-82  picture, 
however,  indicates  a  high  amplitude  wave-like  feature  characterizing  this  flow  both 
at  the  20  m  and  200  m  levels.  It  has  the  appearance  of  being  baroclinically  unstable 
-  especially  at  the  20  m  level.  The  cyclonic  flow  looping  into  the  northeast  of  the 
20  m  map  is  associated  with  the  cold  water  (possibly  upwelled)  in  that  region  shown 
on  the  temperature  map  (Fig.  8a).  It  is  confined  to  the  near-surface  layer  and  can 
be  regarded  as  a  variation  on  the  summertime  picture  (not  shown,  but  stated  by 
Ovchinnikov  (1966)  to  be  basictdly  the  same  as  the  wintertime  circulation). 

On  the  basis  of  the  dynamic  anomaly  map  (Fig.  11a)  the  surface  geostrophic  current 
near  the  mooring  sites  (40°00'N,  006°30'E)  was  southward  at  about  7  cm/s.  This 
geostrophic  flow  itself  does  not  give  rise  to  advective  changes  in  the  mixed-layer 
depth  since  the  flow  will  be  along  isopycnal  surfaces.  Although  other  flow  acting 
near  the  surface,  such  as  Ekman  transport  and  inertial  oscillations,  could  produce 
advective  changes  to  the  mixed  layer,  the  spatial  gradients  near  the  moorings,  and 
for  a  large  distance  surrounding  them,  appear  to  be  very  weak.  Advection  effects  on 
the  mixed-layer  depth  are  therefore  assumed  not  to  be  important  during  MILEX-82. 


3.2.  INTERMEDIATE  LAYER  1150-1000  m) 

Included  less  for  their  direct  relevance  to  the  sea  surface  mixed  layer  but  more 
for  their  general  interest  value  and  for  completeness,  are  south-north  and  west- 
east  transects  of  temperature  (Fig.  12a,b,  respectively)  and  of  salinity  (Fig.  13a,b, 
respectively)  for  the  deeper  150  m  to  1000  m  layer.  These  show  the  presence  of  the 
relatively  salty  (38.5  ppt),  warm  (13.4  °C)  Levantine  Water  between  the  depths  of 
about  300  and  650  m;  a  layer  slightly  deeper  and  thinner  than  that  reported  in  the 
Ligurian  Sea,  for  example,  by  Crepon,  Wald  and  Monget  (1982). 


In  the  north-south  direction  this  Levantine  layer  decreases  in  both  temperature  and 
salinity.  The  temperature  decrease  (Fig.  12a)  is  from  about  13.5  °C  to  13.3  °C,  and 
the  stilinity  decrease  (Fig.  13a)  is  from  38.55  ppt  to  38.50  ppt.  Such  fine  anomalies 
apart,  comparison  of  these  transects  with  portions  of  similar  transects  presented  in 
the  Mediterranean  Atlas  of  Miller,  Tchernia  and  Charnock  (1970)  show  good  general 
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agreement  for  this  layer.  Small  variations  '<re  only  to  be  expected  since  the  so-called 
Levantine- Water  layer  is  highly  variable  throughout  the  Mediterranean. 

Some  weak  temperature  anomalies  are  also  evident  within  the  Levantine  Water  for 
the  west-east  section  (Fig.  12b)  but  far  more  striking  in  this  section  are  the  weirm 
and  cold  temperature  ‘blobs’  sitting  just  above  the  Levantine  layer.  However,  these 
are  largely  density  compensated  by  salinity,  and  any  remaining  density  anomalies 
could  have  no  significant  dynamic  effect. 

■  3.2.1.  Step  structure 

Several  of  the  CTD  profiles  taken  to  the  west  of  the  MILEX-82  zone  displayed  a 
pronounced  step-like  structure  occurring  mostly  below  the  salinity  maximum  and 
within  the  Levantine  layer  as  alternating  relatively  salty  warm  and  fresh  cool  layers. 
As  regards  density,  however,  the  profiles  were  usually  less  step-like.  Two  examples 
of  the  more  pronounced  step-like  structure  are  shown  in  Fig.  14.  In  the  case  of 
CTD  no.  8  an  enlargement  of  the  density  profile  is  shown  with  an  interpretation  in 
terms  of  steps  superimposed  upon  it.  (No  further  discussion  of  these  structtires  is 
presented.) 
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4.  Thermistor  chain  data 


i 


4.1.  ISOTHERM  TIME  SERIES  SECTIONS 


The  thermistor  chain  pairs  (upper  and  lower)  for  each  mooring  were  combined  as 
described  in  Sect.  2.2.  Temperature  contours  were  constructed  by  cubic  spline  in¬ 
terpolation  between  the  22  thermistors  in  the  vertical  and  then  linear  interpolation 
in  time  (i.e.  15  min  depth  values  for  each  isotherm  were  joined  by  straight-line  seg¬ 
ments).  The  resulting  isotherms  are  displayed  as  10-day  sections  in  Fig.  15a, c  for 
mooring  Al,  Fig.  16a, d  for  mooring  B1  and  Fig.  17a,d  for  mooring  Cl.  The  depths 
of  uppermost  and  lowermost  thermistors  are  marked  in  24  h  intervals  on  each  sec¬ 
tion  by  a  star  symbol  and  a  corresponding  instantaneous  temperature  value.  This 
avoided  the  problem  of  labelling  contours  and  also  allowed  a  greater  density  of  them 
to  be  displayed.  Isotherms  are  distinguished  by  simply  counting  in  rounded  half¬ 
degree  intervals  between  printed  temperatures.  As  a  reference  to  mechanical  forcing, 
the  local  wind  speed  is  also  plotted  above  each  section  in  Figs.  15,  16  and  17.  Several 
of  the  thermistors  of  the  lower  thermistor  chain  TC570  malfunctioned  for  much  of 
the  35  days  of  measurement.  These  ‘bad’  data  were  replaced  before  contouring  by 
cubic  spline  interpolations  in  depth  by  utilizing  all  the  data  points  considered  to  be 
‘good’.  Although  great  care  was  taken  in  interpolation  the  spiky  appearance  of  the 
resulting  isotherms  (Fig.  17a,d)  should  be  treated  with  some  caution  -  especially 
since  the  other,  relatively  nearby,  thermistor  chains  A  and  B  do  not  display  the 
same  characteristic. 

A  limitation  of  the  simple  contouring  method  employed  is  its  inability  to  cope  with 
temperature  inversions,  although  these  were  thought  from  preliminary  examination 
to  be  unlikely,  especially  in  the  upper  thermocline.  The  method  was  to  search  for 
isotherms  begitming  at  the  top  suid  terminate  search  at  the  first  occurrence  of  the 
required  temperature  value.  If  temperature  inversions  were  present  at  mooring  C 
they  might  also  account  for  the  confusing  shape  of  the  isotherms  there.  It  is  unlikely, 
however,  that  inversions  would  occur  only  at  C  and  not  at  A  and  B 
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4.2.  INERTIAL  OSCILLATIONS 

We  see  irom  the  isotherm  displacement  in  Figs.  15  to  17  that  the  thermocline  is 
characterized  by  strong  inertial  waves.  The  wave  amplitudes  were  typically  about 
3  m  or  4  m  throughout  most  of  the  35  days  of  measurements.  They  were  greatest 
during  stormy  conditions  following  rapid  increases  in  the  wind  speed  such  Ets  oc¬ 
curred  around  days  279  and  280  and  they  appear  to  have  been  greater  at  mooring  A 
(Fig.  15b)  than  at  the  other  two  moorings  (cf.  Figs.  16b  juid  17b).  All  the  plots  show 
a  tendency  for  the  amplitude  to  increase  with  depth;  a  good  example  being  during 
the  relatively  calm  weather,  such  as  days  270  to  280  (Fig.  15a),  when  inertial  waves 
can  hardly  be  seen  in  the  highly  stratified  water  of  the  upper  thermocline  but  are 
fairly  pronoimced  lower  down  (amplitude  of  up  to  5  m).  The  same  effect  occurred 
during  the  stormier  conditions  around  days  280  to  284  (Fig.  15b)  with  amplitude 
increasing  through  the  thermocline  from  about  6  to  edmost  9  m.  The  total  energy 
of  the  vertical  oscillations,  however,  may  not  be  increasing  with  depth  in  the  same 
way  since  the  buoyancy  restoring  forces  due  to  stratification  are  far  higher  in  the 
upper  than  in  the  lower  thermocline. 

There  appears  to  be  a  significant  delay  between  increases  in  wind  speed  and  the 
onset  of  vertical  oscillations.  A  good  example  occurred  between  days  278  and  281  at 
mooring  A  (Fig.  15a,b).  The  wind  increased  fairly  rapidly,  starting  about  6  h  before 
day  279  at  1  m/s  and  reaching  16  m/s  9  h  later.  No  corresponding  response  can 
be  observed  in  the  thermocline  at  this  stage.  Over  the  next  18  h  the  wind  dropped 
back  to  2  m/s  and  increased  again  to  13  m/s.  Only  then,  some  24  h  after  the  onset 
of  these  strong  wind  changes  did  the  isotherms  begin  to  respond.  The  oscillations 
began  with  a  strong  downward  displacement  of  the  isotherms  and  then  continued 
to  oscillate  at  near  inertial  frequency  for  about  6  days;  by  which  time  they  had 
attenuated  considerably. 


4.3.  GENERATION  OF  VERTICAL  OSCILLATIONS 

Why  should  there  be  such  a  delayed  response  of  the  thermocline  to  the  winds?  The 
answer  will  depend  on  which  generation  mechanism  causes  the  vertical  oscillations. 


■  4.3.1.  Horizontal  advection 

The  possibility  of  horizontal  advection  of  sloping  isotherms  by  currents  can  be  dis¬ 
counted  as  a  mechanism  since  there  is  no  evidence  of  even  weak  thermal  fronts 
during  MILEX-82  (see  Sect.  3.1).  Indeed,  were  this  to  be  the  mechanism,  we  note 
that  taking  inertial  circles  to  be  of  order  4  km  diameter  at  the  thermocline  depth  (a 
typical  value  taken  from  progressive  vector  diagrams  not  shown  here)  and  tempera¬ 
ture  changes  of  up  to  4  °C  per  inertial  cycle  (Fig.  15b),  a  front  as  strong  as  1  °C/km 
would  be  necessary  to  fully  account  for  the  observed  vertical  oscillations  through 
advection.  In  any  case,  since  little  delay  occurs  between  sudden  wind  changes  and 
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the  horizontal  oscillations  (see  Fig.  23  of  Sect.  5.1)  such  advective  effects  could  not 
explain  the  delay  in  vertical  oscillations. 

The  importance  of  advection  should  not  be  de-emphasized  too  strongly,  however, 
since  there  is  evidence  from  CTD  data  for  spatial  intermittance  of  temperature 
emd  salinity  structures  in  the  lower  thermocline  similar  to  those  in  the  Levantine 
Water  shown  in  Fig.  14  of  Sect.  3.  If  such  structures  were  present  in  the  vicinity  of 
mooring  C,  for  example,  their  advection  back  and  forth  by  the  inertial  currents  might 
explain  the  abrupt  changes  and  ‘top-hat’  shapes  observed  in  the  isotherms  there 
(Fig.  17a,d  -  see  also  cautionary  remarks  concerning  accuracy  of  these  isotherms  in 
Sect.  4.1).  The  shape  of  the  vertical  oscillations,  at  least,  will  certainly  be  modified 
by  advective  effects,  although  they  do  not  constitute  vertical  motion. 

There  are  essentially  two  theories  on  the  direct  forcing  of  vertical  inertial  oscilla¬ 
tions  and  the  consequent  downward  propagation  of  inertio-gravity  waves.  Both  of 
these  could  account  for  the  observed  delay  in  thermocline  response  to  the  winds. 
The  first,  due  to  KroU  (1975),  requires  a  mean  Ekman  mass  flux  divergence  and 
hence  spatial  variations  in  the  wind  stress.  The  second,  due  to  KroU  (1982)  and 
KroU  (1984),  involves  a  parametric  instabiUty  of  the  slab-Uke  mixed-layer  flow  and 
does  not  depend  on  spatial  variations  in  wind  stress. 


■  4,3.2.  Spatially  varying  wind  stress 

To  iUustrate  in  a  simple  way  how  spatial  variations  in  wind  stress  can  induce  vertical 
osciUations  we  write  down  the  Unearized,  bulk  (verticaUy  integrated)  momentum 
balance  for  the  mixed  layer  reduced  to  the  form 

^ -(-/ic  X  U-1  ^To  =  0,  (1) 

where  U  tmd  V  are  the  bulk  horizontal  velocity  components  in  an  orthogonal  coor¬ 
dinate  system  {x,!/,z)  with  z  directed  (-l-w)  upward,  /  is  the  CorioUs  parameter,  k 
is  the  unit  vertical  vector,  and  t*  and  Vy  are  the  surface  wind  stress  components.  A 
solution  to  Eq.  (1)  for  a  suddeiUy  imposed  wind  stress  describes  the  famiUar  inertial 
osciUations  and  Ekman  mass  transport  for  the  mixed  layer  (PoUmU  et  al.,  1973). 
Taking  the  divergence  and  curl  of  Eq.  (1)  and  using  the  continuity  relation 


dU  dV 
ax  ay 

where  z  =  -fi  is  the  mixed-layer  base,  yields 


(2) 


'*  -  /  curl  U  -1  -  div  To  =  0 
at  p 

d  1 

— (curlU)  -I  +  -  curl To  =  0. 

at  p 


(3) 
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A  solution  to  Eq.  (3)  for  a  suddenly  imposed  wind  stress  which  then  remains  constant 
in  time  but  varies  spatially  is  given  by 

u)_a  =  — divTosin/(  +  curlro(l  -  cos  ft), 

IP 


curlU  =  —  divro(l  -  cos  ft)  -  curl  To  sin /t. 


(4) 


We  see  from  Eq.  (4)  that  vertical  oscillations  can  arise  due  to  either  a  purely  diver¬ 
gent  or  a  purely  rotational  wind  stress.  In  addition,  wind-stress  divergence  imparts 
mean  vorticity  to  the  mixed  layer  whereas  wind-stress  curl  causes  the  familiar  Ek- 
man  suction  velocity  at  the  mixed-layer  base.  The  latter  was  invoked  by  Kroll  (1975) 
to  force  inertial  oscillations  in  an  inviscid  stratified  interior  which  he  then  examined 
using  ray  theory.  Rubenstein  (1983)  used  a  similar  approach  but  did  not  decouple 
the  interior  from  the  mixed  layer,  thus  allowing  an  interaction  between  the  two. 
Also,  using  a  numerical  model,  Rubenstein  treated  all  frequency  components  of  the 
forcing  simultaneously  rather  than  specifying  a  single  frequency  component.  Ear¬ 
lier  models  by  Pollard  (1970)  and  Pollard  and  Millard  (1970)  used  the  same  forcing 
mechanism  (wind-stress  curl)  emd  examined  the  inertial  oscillations  as  a  function 
of  depth  by  means  of  normal  mode  expansions  in  the  vertical.  Krauss  (1979)  ad¬ 
dressed  the  question  of  inertial  oscillations  (horizontal  and  vertical)  in  a  channel  of 
finite  width  and  depth  but  infinite  in  length.  Krauss’  model  was  applied  first  to 
one  homogeneous  layer  and  then  to  two  homogeneous  layers  of  different  densities. 
Two  types  of  forcing  were  applied  in  each  case;  firstly,  using  a  purely  divergent  wind 
stress  and  secondly,  using  a  purely  rotational  wind  stress. 

The  above-mentioned  models  differ  from  each  other  in  their  various  refinements,  in 
their  methods  of  solution  and  in  their  emphasis  on  the  controlling  parameters.  The 
basic  forcing  mechanism,  however,  is  common  to  all.  The  inertial  currents  must 
be  generated  by  a  storm  which  is  botmded  horizontally  and  which  therefore  has 
associated  winds  which  vary  spatieJly.  This  is  doubly  significant  since  not  only  will 
near-inertial  vertical  oscillations  (Eq.  (4))  develop  at  the  pycnocline  and  mixed-layer 
interface,  thereby  constituting  an  energy  source  for  inertial  waves  in  the  pycnocline, 
but  it  will  also  establish  a  spatial  scale  defining  the  horizontal  wavenumber  k  in  the 
dispersion  relation  for  internal  waves: 

=  f^  +  k‘ I  ,  (5) 

where  N  is  the  Brunt-Vaisala  frequency  and  m  the  vertical  wavenumber.  The  ver¬ 
tical  group  velocity  associated  with  Eq.  (5)  is 

c,  =  -k^  Iwm^.  (6) 

It  was  shown  quite  rigorously  in  Kroll  (1975)  that  the  forcing  frequency  at  the  mixed- 
layer  base  will  differ  from  the  local  value  of  /  by  a  small  positive  amount  if  pressure 


SaCLANTCEN  SM-213 


tenns  are  included  in  the  Ekman  layer  equations;  in  our  case  Eq.  (1).  Kroll  also 
points  out  that  the  non-linear  terms  excluded  from  the  same  equations  may  also 
have  a  modifying  influence  on  the  inertial  frequency.  It  seems  reiisonable  to  suppose 
the  same  will  be  true  of  horizontal  diffusion. 

Since  the  forcing  frequency  w  is  determined  in  the  surface  layer  and  must  be  greater 
than  /,  and  since  k  is  also  determined  by  the  spatial  scale  of  the  storm,  Eq.  (5) 
shows  that  m  is  finite  and  the  vertical  group  velocity  (Eq.  (6))  is  non-zero  and 
positive  downward.  Near-inertied  energy  thus  propagates  downward  at  small  angles 
to  the  horizontal  and  can,  in  the  relatively  low  levels  of  turbulent  dissipation  in 
and  below  the  pycnocline,  propagate  over  long  distances.  If  this  mechanism,  one  of 
strong  spatial  variations  in  the  wind-stress  forcing,  is  in  fact  operating  in  or  near 
the  MILEX-82  area,  one  would  expect  the  forcing  to  take  place  in  the  areas  of 
strongest  wind-stress  curl  several  hundreds  of  kilometres  to  the  north.  This  area  is 
known  (May,  1982)  to  have  an  exceptionally  high  wind-stress  curl  both  in  monthly 
and  seasonal  means  due  to  persistence  of  the  Mistral  wind  with  its  narrow  jet¬ 
like  profile.  In  support,  meastued  values  of  the  Ekman  suction  velocites  have  been 
reported  (Voorhis  and  Webb,  1970)  during  the  MEDOC  cruise  of  up  to  2  cm/s. 
Such  events  would  be  accompsuiied  by  extremely  energetic  near-inertial  oscillations 
in  the  vertical. 

A  possible  explanation  for  the  delayed  response  of  the  thermocline  to  the  wind, 
therefore,  might  be  that  the  response  is  not  to  the  local  wind  but  rather  to  the  wind 
250  km  to  the  North.  The  inertial  waves  would  have  to  be  travelling  southward  with 
group  velocity  about  250  km  in  18  h.  The  current  meters,  however,  in  contrast  to 
their  thermistor  chains,  do  not  show  this  same  delayed  response  to  the  wind  (see 
Fig.  23  of  Sect.  5)  and  so  do  not  support  this  distant  source  hypothesis.  On  the 
other  hand,  as  pointed  out  by  Schott  (1971),  the  most  energetic  waves  for  the  hori¬ 
zontal  motions  (detectable  by  current  meters)  may  be  separated  in  the  wavenumber 
spectrum  from  the  most  energetic  waves  for  the  vertical  motions  (detectable  by  ther¬ 
mistor  chains);  the  vertical  motions  tending  to  associate  with  higher  wavenumber. 
In  this  case  one  might  expect  a  different  group  velocity  and  arrival  time  for  the 
respective  wave  trains. 


■  4.3.3.  Instability  mechanisms 

More  recent  theories  for  downward  propagating  inertial  waves  involve  instability 
mechanisms  due  to  Kroll  (1982)  and  Kroll  (1984).  In  Kroll  (1982)  a  uniform  slab 
flow  of  the  mixed  layer  is  assumed,  in  which  there  are  no  inertial  oscillations.  The 
linearized  perturbation  equations  for  the  upper  slab  layer  and  the  lower  stratified 
layer  are  shown  to  be  unstable  for  certain  combinations  of  the  governing  parame¬ 
ters.  The  length  scale  of  the  instability  is  of  the  order  10  km  and  is  thus  supportive 
of  observations  of  coherence  scales  for  inertial  motions  below  the  ntixed  layer  (Pol¬ 
lard,  1980;  Sandford,  1975).  An  interesting  feature  of  the  instability  is  that  it  does 
not  depend  on  spatial  variations  in  wind  stress. 
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The  second,  related,  instability  mechanism  (KroU,  1984)  is  both  more  realistic  (since 
it  includes  inertial  oscillations  in  the  ‘mean’  slab-flow)  and,  as  it  turns  out,  more 
unstable.  The  linearized  perturbation  equations  are  essentially  those  of  KroU  (1982) 
but  with  coefficients  that  include  time  periodic  components  (i.e.  components  repre¬ 
senting  the  inertial  currents  in  the  mixed  layer).  Instability  au'ising  in  such  equations 
is  known  as  parametric  instability  and  the  solutions,  especiaUy  in  this  ceise,  are  rather 
complicated.  Nevertheless,  some  interesting  results  which  are  relevant  to  the  inter¬ 
pretation  of  the  MILEX-82  data  emerge.  In  a  comparison  of  theory  with  data  KroU 
chose  parameters  observed  at  site  D  presented  in  PoUard  (1980).  He  found  that 
instabiUty  was  possible  for  a  small  region  of  the  pareuneter  space  provided  the  hori¬ 
zontal  length  scale  L  =  V^h/ /  I®*®  about  2.5  km.  Here  i/h  is  a  coefficient 

for  horizontal  turbulent  diffusion  in  the  mixed  layer.  To  ensure  a  good  correspon¬ 
dence  with  the  measured  mean  flow,  however,  a  rather  smaU  choice  of  £  must  be 
made  which  leads  to  a  horizontal  wavelength  for  the  near  inertial  osciUations  as  low 
as  1.2  km.  Although  this  seems  unrealisticaUy  short  for  a  wavelength,  it  represents 
a  lower  botmd,  and  the  estimate  in  any  case  appears  quite  sensitive  to  the  choice  of 
some  of  the  variable  parameters.  Since  there  is  a  reasonable  similarity  between  the 
PoUard  (1980)  data,  referred  to  by  KroU,  and  that  of  MILEX-82,  we  might  expect 
the  instability  to  be  operating  in  the  MILEX-82  area.  To  be  more  precise  one  would 
have  to  make  fresh  calculations  using  KroU’s  model. 

Another  set  of  parameters  presented  in  KroU  (1984),  as  realistic,  and  in  our  case 
bearing  some  similarity  to  MILEX-82  data,  gave  an  e-folding  growth  time  for  the 
most  unstable  solution  of  1.3  days.  This  growth  time  could  offer  an  alternative 
explanation  of  the  delayed  thermocline  response  to  the  rapid  changes  in  wind  stress. 

In  summary,  it  is  Ukely  that  both  the  spatiidly  varying  wind  stress  and  the  paramet¬ 
ric  instability  mechanisms  could  explain  the  presence  of  vertical  osciUations  in  the 
MILEX-82  area.  Furthermore,  both  mechanisms  can  explain  the  delayed  response 
of  the  thermocUne  to  the  wind.  To  determine  which  is  the  more  Ukely  mechanism 
critical  tests  would  be  requited  applying,  for  example,  KroU’s  (1984)  model  to  the 
MILEX-82  parameters  on  the  one  hand,  and  then  on  the  other  hand  applying  either 
the  model  due  to  KroU  (1976)  or  Rubenstein  (1982)  to  the  hypothesized  generation 
zone  to  the  North  and  studying  the  effect  in  the  MILEX-82  area. 


4.4.  SUPER-INEPTIAL  FREQUENCIES 


It  is  quite  evident  from  the  isotherm  plots  of  Figs.  15  to  17  that  vertical  osciUa¬ 
tions  are  also  occiuring  at  approximately  twice  the  local  inertial  frequency.  A  good 
example  is  again  the  period  covering  days  280  to  286  in  Fig.  15b.  By  day  282  a 
very  pronounced  double  peak  is  occurring  on  the  inertial  waves  and  as  time  passes 
the  inertial  motions  attenuate  and  become  confused  with  the  super-inertial  motions 
(w  Si  2f)  which  themselves  remain  detectable  until  the  end  of  the  record  (Fig.  16c). 
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Similar  observations  of  super-inertial  motions  have  been  reported  by  Titov  (1973), 
Krauss  (1981)  and  Pinkel  (1983).  Although  not  specifically  mentioned,  other  exa¬ 
mples  can  be  found  in  the  temperature  variance  spectra  of,  for  example,  Weller 
euid  Halpern  (1983)  and  Moen  (1984). 

Relative  to  the  inertial  waves,  these  higher  frequency  components  are  associated  with 
high  downward  vertical  group  velocities  (see  Eqs.  (5)  and  (6)),  thus  making  them 
efficient  in  extracting  energy  from  the  mixed-layer  flow.  This  has  been  pointed 
out  by  several  investigators,  for  example  Krauss  (1981),  Rubenstein  (1983)  and 
KroU  (1984), 


4.5.  TEMPERATURE  SPECTRA 

To  better  illustrate  the  presence  of  super-inertial  verticeJ  oscillations,  spectra  of 
temperature  variance  are  presented  in  Figs.  18  to  20  and  spectra  of  isotherm  dis¬ 
placement  in  Fig.  21.  In  these  (particularly  in  Fig.  21)  significant  pe2>ks  can  be  seen 
at  roughly  integral  multiples  of  the  inertial  frequency;  in  some  cases  up  to  a  factor 
of  four.  Rubenstein  (1983)  has  pointed  out  that  such  super-inertial  components 
are  probably  the  result  of  non-linear  interaction  between  near-inertial  waves  and 
represent  a  sink  of  energy  for  the  inertial  motions.  In  Kroll  (1984)  super-inertial 
waves  are  a  prediction  of  the  parametric  instability  mechanism  and  can  be  taken, 
according  to  Kroll,  as  evidence  that  the  mechanism  is  operating.  The  frequencies 
in  Kroll’s  case  should  occur  not  at  n/  but  rather  at  u)  4-  n/,  where  n  is  an  integer 
and  u)  is  the  near  inertial  frequency  of  the  unstable  oscillations  in  the  mixed-layer 
flow.  Unfortunately,  due  averaging  over  five  adjacent  frequency  bands,  the  resolu¬ 
tion  of  the  spectra  in  Figs.  18  to  20  is  not  good  enough  to  distinguish  between  n/ 
suid  w  +  n/  for  the  super-inertial  frequencies.  For  such  a  fine  distinction  the  records 
ought  to  be  re-anedysed  using  complex  demodulation. 


4.6.  PHASE  RESOLUTION 

Before  continuing  to  discuss  current-meter  data,  it  is  interesting  to  examine  the 
isotherm  plots  for  phase  changes  with  depth.  In  the  case  of  shallow  seas,  for  example, 
it  has  been  both  observed  and  predicted  (Krauss,  1979)  that  a  two-layer  structure 
to  the  inertial  motions  will  develop  with  an  upper  layer  180°  out  of  phase  with  a 
lower  layer.  It  i.s  not  clear  that  the  same  should  be  true  of  the  vertical  oscillations 
and  especisdly  in  much  deeper  water  such  as  in  the  MILEX-82  area.  The  clearest 
example  of  phase  changes  that  can  be  seen  directly  in  the  isotherm  plots  is  for 
days  280  to  286  at  mooring  A.  The  phase  shifts,  illustrated  in  Fig.  22,  indicate  that 
the  phase  does  change  over  some  region  of  the  thermocline  between  40  and  60  m 
but  the  change  is  not  constant  in  time  and  does  not  reach  the  value  180°  (i.e.  9.3  h 
at  this  latitude).  The  upper  and  lower  levels  of  the  thermocline  appear  to  be  in 
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phase  at  the  onset  of  the  vertical  oscillations  on  day  280.  Curiously,  the  upper  level 
then  begins  to  lag  the  lower,  first  by  a  few  hours,  and  reaching  a  maximum  of  7  h 
by  day  285.  Although  the  vertical  group  velocity  is  downward,  therefore,  as  defined 
by  Eq.  (6),  Fig.  22  indicates  an  upward  directed  vertical  phase  speed  (Cp  =  ui/m) 
consistent  with  observations  by  Pollard  (1980).  Whether  or  not  the  phase  changes  of 
Fig.  22  are  typical  of  all  the  thermistor  chain  data  is  difficult  to  tell  without  further 
analysis  such  as  complex  demodulation  mentioned  earlier. 
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Of  particular  interest  in  Fig.  25  are  phase  angles  between  the  40  and  60  m  depths 
and  the  60  and  80  m  depths,  since  these  three  current  meters  were  on  the  same 
mooring,  until  day  295  when  the  mixed-layer  deepened  to  beyond  40  m  (Fig.  31), 
they  were  situated  well  within  the  pycnocline.  On  day  270  the  40  and  60  m  inertial 
currents  were  in  phase,  but  by  day  275  they  were  out  of  phase  by  about  245°,  where 
they  remeuned  for  about  4  to  5  days.  They  came  back  into  phase  by  day  282.5,  then 
went  about  180°  out  again  during  the  period  spanning  days  285  to  287.5,  back  into 
phase  by  day  290,  90°  out  by  day  295.5  and  fmally  remaining  within  45°  for  the 
last  10  days.  The  recordings  for  6f  and  80  m  remain  closer  in  phase  apart  from  a 
brief  period  around  day  287.5  when  they  become  almost  180°  out  of  phase.  These 
periods  of  antiphase  between  closely  spaced  current  meters  would  tend  to  raise  the 
level  of  turbulent  diffusion  through  shear  production.  The  levels  of  shear,  together 
with  buoyancy  frequency  and  Richardson  number  are  presented  in  Sect,  6. 

In  other  examples  of  phase  variations  which  can  be  followed  in  Fig.  25  it  would  not 
be  possible,  due  to  horizontal  separation  of  current  meters,  to  attribute  phase  dif¬ 
ferences  oidy  to  vertical  dispersion.  Nevertheless,  the  high  variabiUty  in  phase  over 
such  small  separations  is  surprising  and  suggests  that  coherence  between  records 
estimated  over  the  entire  measuring  period  would  be  low.  Although  it  will  not  be 
presented  here,  a  more  thorough  treatment  would  require  estimating  phase  as  a  func¬ 
tion  of  time  by,  for  example,  the  method  of  complex  demodulation  (Perkins,  1970; 
Pollard  and  Tarbell,  1975).  This  would,  in  addition,  provide  information  on  the 
small  shifts  in  frequency  from  the  local  inertial  value. 


S.3.  FREQUENCY  DISPERSION 


The  variations  in  phase  described  above  (Sect.  5.2)  might  be  due  in  part  to  a  vari¬ 
ation  in  frequency  of  the  near-inertial  motions  with  depth.  Although  complex  de¬ 
modulation  was  not  performed  on  the  MILEX-82  data,  some  idea  of  the  variations  in 
near  inertial  frequency  can  be  extracted  from  the  stick  plots  of  Fig.  23.  The  records 
for  each  current  meter  were  divided  into  5-day  intervals  and  the  mean  inertial  period 
for  each  interval  was  estimated  by  measuring  the  time  required  for  5,  6  or  7  complete 
oscillations.  Cases  that  were  not  clear,  or  that  contained  sudden  phase  changes  were 
not  included.  End  points  were  determined  by  requiring  end  vectors  to  be  tdigned,  if 
necessary,  by  an  approximate  interpolation  between  adjacent  vectors.  Results  of  this 
analysis  are  presented  in  Table  1  as  a  percentage  above  the  local  inertial  frequency 
(/  =  2)r/18.66  h-'-). 

We  see  that  frequency  is  indeed  dispersed  in  some  way  both  with  depth  and  time  - 
the  extremes  ranging  from  -2%  to  -hl5%.  Of  particular  interest  are  the  differences 
between  the  40  and  60  m  depths  demonstrated  earlier  to  be  variously  in  and  out 
of  phase.  It  would  appear  from  this  somewhat  crude  analysis  that  the  frequency  at 
40  m  is  decreasing  with  time  from  about  1.15/  on  day  275  to  about  /  by  day  300.  At 
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Table  1 

Mean  inertial  periods  for  selected  5-day  intervals 


D*‘pth 

(m) 

Time  interval  (Julian  days) 

Mean 

inertial 

period 

270/ 

275 

275/ 

280 

280/ 

285 

285/ 

290 

290/ 

295 

295/ 

300 

300/ 

305 

to 

0 

0 

2 

0.7 

20 

0 

_  1 

2 

0.3 

30 

7 

4 

4 

0 

4.3 

40 

10 

15 

4 

4 

2 

0 

1 

5.1 

60 

5 

5 

7 

5 

15 

2 

6.5 

80 

5 

-  1 

5 

7 

0 

4 

2 

3.1 

120 

2 

2 

13 

8 

^  2 

4.6 

160 

4 

4 

4 

8 

15 

4 

6.5 

60  m,  on  the  other  hand,  the  frequency  increases  over  the  same  period  from  about 
1.05/  to  1.15/.  Looking  at  the  mean  values,  there  is  a  tendency  for  the  frequency 
to  increase  with  depth,  or  at  least  for  the  mixed-layer  oscillations  to  be  significantly 
lower  in  freqeuncy  (<  1.01/)  than  those  of  the  pycnocline  (1.05/). 

This  sort  of  dispersion  can  be  explained  by  the  dispersion  relation  (Eq.  4.5)  and 
vertical  group  velocity  expression  (Eq.  4.6)  for  inertial-internal  waves.  As  pointed 
out  by  several  authors  (e.g.  Kroll,  1975;  Krauss,  1979;  Rubenstein,  1983)  sudden 
temporal  changes  in  wind  speed  can  excite  a  spectrum  of  near-inertial  waves  in  the 
upper  ocean,  and  for  a  fixed  stratification  N  those  wave  packets  with  lowest  vertical 
wavenumber  m  will  have  the  highest  frequency  w  .  Also,  these  same  wave  packets  will 
have  the  highest  downward  group  velocities,  causing  the  higher- frequency  inertial 
energy  to  disappear  first  from  the  upper-most  layer.  Of  course  the  mixed  layer  is  not 
stratified  and  so  for  finite  vedues  of  where  k  is  the  horizontal  wavenumber, 

Eq.  (4.5)  shows  that  the  frequency  will  tend  towards  /. 

There  are  other  dynamics-related  reasons  why  the  inertial  frequency  just  below  the 
mixed  layer  will  be  greater  than  /.  These  were  discussed  in  Sect.  4.3  in  cotmection 
with  the  vertical  oscillations.  It  was  pointed  out  that  the  oscillations  can  only  enter 
the  pycnocline  below  the  mixed  layer  through  either  an  Ekman  suction  mechanism 
(KroU,  1975)  or  through  instabiUty  mechanisms  of  the  mixed-layer  flow  (Kroil,  1983; 
Kroll,  1984).  It  was  also  argued  that  if  the  deeper  oscillations  were  due  to  the  for¬ 
mer  mechanisms  they  were  probably  not  generated  locally,  but  rather  to  the  north 
in  regions  of  high  spatial  variations  in  wind  stress.  Support  for  the  instability  mech¬ 
anisms  for  the  generation  of  sub-mixed-layer  oscillations  comes  with  the  observation 
of  energy  peaks  in  the  velocity  spectra  (Fig.  26)  at  higher  harmonics  of  the  in- 
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ertial  frequency  (KroU,  1984).  These  were  also  observed  in  the  temperature  and 
isotherm  spectra  presented  in  Sect.  4.5.  Although  such  super-inertial  oscillations 
are  a  prediction  of  Kroll’s  instability  theory,  inertied  oscillations  are,  in  any  case, 
weakly  non-linear  and  so  ^ls  pointed  out  by  Rubenstein  (1984),  a  non-linear  term 
Uke  u  dv/dx  can  easily  reach  10%  of  linear  terms  in  the  momentum  balance.  The 
effect  of  non-linearities  is  to  generate  higher  harmonics  of  the  inertial  frequency. 
Another  strong  non-linearity,  of  course,  occurs  in  the  transition  zone  between  mixed 
layer  and  pycnocline.  Here  one  can  compare  the  i-momentum  balance  term  w  dujdz 
with  fv,  assuming  the  pycnocUne  to  be  initally  quiescent;  if  not,  the  inertial  motions 
can  he  out  of  phase  with  the  mixed  layer  as  is  often  the  case  (Krauss,  1979;  also 
recall  the  phase  lag  between  20  and  30  m  noted  earlier  in  Sect.  5.2)  causing  even 
higher  values  of  the  shear  dujdz.  Estimates  of  the  shear  in  the  20  to  30  m  layer  (see 
Fig.  27  of  Sect.  6)  show  sustained  values  of  0.05  s"* .  For  the  same  period  (days  280 
to  285),  the  isotherm  plots  of  Fig.  15  show  a  crest-to- trough  range  near  the  interface 
of  4  m/inertial  cycle,  which  converts  to  a  vertical  velocity  amplitude  of  0.C2  cm/s. 
Taking  a  t/-velocity  of  15  cm/s,  which  is  about  the  amplitude  deduced  by  linear 
interpolation  between  20  and  30  m  or  between  30  and  40  m  (Fig.  23),  gives  a  value 
of  1.5  X  10"^  cm/s^  for  fv  compared  to  1.0  x  10“^  cm/s^  for  w  dujdz.  Non-Unear 
terms,  therefore,  are  of  the  same  order  as  Unear  terms  and  it  is  hardly  surprising 
that  super-inertial  osciUations  are  observed. 

The  discussion  to  this  point  has  been  on  near-inertial  frequencies,  including  higher 
harmonics,  which  occur  below  the  mixed  layer.  It  is  important  to  note  that  the 
essentially  horizontal  oscillations  which  characterize  the  mixed  layer  can  also  differ  by 
small  amounts  in  frequency  from  the  local  value  of  /.  There  appear  to  be  three  main 
reasons  why  this  is  so.  The  first  reason  has  to  do  with  the  effect  of  small  pressure 
gradient  forces  (KroU,  1975)  usuaUy  neglected  in  boundary  layer  analysis  of  inertial 
waves.  KroU  showed  analyticaUy  how  pressure  effects  in  the  Unearized  boundary 
layer  equations  can  cause  a  sUght  increase  in  the  inertial  frequency  and  then  he 
assumed  the  result  to  hold  in  the  non-linear  case.  On  the  basis  of  Kroll’s  analysis, 
and  as  mentioned  earUer  in  Sect.  4.3,  one  might  assume  further  that  other  effects,  also 
normaUy  neglected,  such  as  Unear  and  non-linear  advection  and  horizontal  diffusion, 
would  also  cause  smaU  changes  in  the  inertial  frequency. 

The  second  reason  for  observed  changes  in  the  inertial  frequency  could  be  due  to  a 
doppler  shift.  Schott  (1971)  has  pointed  out  that  for  inertial  waves  with  a  50  km 
wavelength  such  as  he  observed  in  the  northern  North  Sea,  a  1%  frequ  ncy  variation 
near  the  inertial  frequency  oiUy  needs  a  mean  current  of  1  cm / s  along  the  propagation 
path  of  the  wave. 

Finally,  the  third  reason  involves  spatial  variations  in  the  underlying  quasi-geostro- 
phic  flow  field  (Fomin,  1973;  Weller,  1982).  Contributions  by  other  authors  to  this 
topic  Me  reviewed  in  WeUer’s  paper.  In  simple  terms,  water  particles  in  the  mixed 
layer  accelerated  by  changes  in  wind  stress  attempt  to  move  in  straight  lines  confined 
to  the  earth’s  surface  by  gravity.  However,  since  the  ‘background’  is  rotating  both 
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due  to  et^rth  rotation  and  pre-existing  vorticity,  the  particles  will  appear  to  travel 
in  elliptical  orbits.  The  frequency  of  these  motions  will,  therefore,  differ  from  /,  the 
local  inertial  frequency  fixed  by  latitude,  according  to  the  sign  and  strength  of  the 
background  vorticity. 

With  reference  to  the  dynamic  anomcdy  plots  of  Fig.  11  and  the  stick  diagrams  of 
Fig.  23,  both  the  doppler  rind  vorticity  effects  on  /  ought  to  he  minimal  in  the 
MILEX-82  region.  This  might  explain  the  relatively  smcdl  deviations  to  /  presented 
in  the  mixed  layer  portion  of  Table  1 . 


5.4.  velocity  spectra 

Rotary  spectra  for  the  upper  6  current  meters  are  presented  in  the  composite  plot 
of  Fig.  26.  Those  for  the  NBA  current  meters  at  120  m  and  160  m  have  been 
omitted  since  there  is  evidence  that  these  were  severly  ‘pumped’  by  mooring  line 
oscillations  in  the  vertical  as  discussed  earlier  in  Sect.  2.2.  The  method  of  averaging 
over  adjacent  frequency  hands,  as  a  method  of  increasing  confidence  in  the  spectra, 
has  decreased  resolution  to  the  extent  that  small  deviations  from  /  in  the  local 
inertial  frequency  caimot  be  resolved. 

The  spectra  of  Fig.  26  confirm  the  predominance  observed  earlier  in  the  stick  plots 
(Fig.  23)  of  the  highly  clockwise  rotational  inertial  oscillations.  Not  only  is  there 
a  peak  at  the  inertial  frequency  but  clockwise  energy  appears  to  saturate  a  large 
band  of  super-inertial  frequencies  down  to  about  4  h  periods.  Within  this  band, 
as  with  the  temperature  variance  spectra  of  Figs.  18  to  20,  smaller  peaks  can  be 
observed  at  higher  harmonics  of  the  near-inertial  frequency.  In  the  mixed-layer 
spectra  ( 10  and  20  m)  these  peaks  must  be  due  to  non-linear  effects  discussed  above 
(Sect.  5.3),  whereas  in  the  deeper  layers  (particularly  60  and  80  m),  where  arguably 
the  clockwise  band  has  relatively  more  energy,  they  could  be  due  to  both  non¬ 
linearity  and  parametric  instability  mechanisms  (Kroll,  1982  and  1984)  discussed  in 
Sects.  4.3  and  5.3.  The  anticlockwise  spectrum  at  40  m  shows  quite  pronoimced 
peaks  at  the  inertial  and  double  inertial  frequencies,  implying  that  the  motion  is 
more  elliptic  at  40  m  than  for  the  other  depths . 
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6.  other  parameters 

This  section  contains  time  series  plots  of  various  measured  or  estimated  parameters 
which  are  directly  relevant  to  mixed-layer  deepening  and  which  will  be  of  interest  in 
any  comparison  of  the  data  with  mixed-layer  models  for  the  MILEX-82  area. 

Three  sets  of  parameters  are  presented.  In  the  first  set,  bulk  values  of  sheeir  5, 
Brunt- Vaistila  frequency  AT,  and  Richardson  number  Ri  were  calculated  according 
to  the  formulae 


5*  = 


(Ul  -  -1  (t),  -  Vjf 

(^1 

(T, -Tj) 


= -go.  .  , 
Ri  =  JVV5^ 


(7) 


where  the  subscripts  1  and  2  refer  to  the  upper  and  lower  current  meter  levels 
respectively,  the  temperatures  Ti  and  Tt  were  estimated  by  linear  interpolation  on 
the  thermistor  chain  data,  and  a  is  the  coefficient  of  thermal  expansion  for  sea  water. 
Errors  in  expression  (Eq.  (7))  occur  due  to  neglecting  effects  of  horizontal  separation 
of  ciu-rent  meters,  neglecting  differences  in  velocity  due  to  type  of  current  meter  and 
in  neglecting  salinity  effects  on  density.  The  latter  error  is  of  the  order  10%  in  the 
MILEX-82  area.  Plots  of  S,  N  and  Ri  are  displayed  in  Figs.  27  to  30.  Values  for  the 
intervals  80  to  120  m  and  120  to  160  m  are  not  shown  since  the  two  NBA  current 
meters  (at  120  and  160  m),  as  discussed  earlier,  produced  unreliable  results. 


The  second  set  of  parameters  (Fig.  31a,d)  are;  mixed-layer  depth,  maximum  Brunt- 
Vaisala  frequency,  change  in  heat  content  and  change  in  potential  energy.  To  calcu¬ 
late  these,  cubic  spline  functions  were  fitted  to  the  thermistor  chain  data  and  then, 
using  the  appropriate  analytic  formula  for  each  spline  section,  the  depth  of  minimum 
second  derivative,  d^Tfdz^  was  defined  as  the  mixed-layer  depth,  the  maximum  first 
derivative,  dT/dz  was  used  to  calculate  N,  and  the  integrals 

AH{t)  =  pC^f  {T(t}-T(0))dz, 

} -D 

^PE(t)  =  -gpal  {T{t)-T(0))zdz  (8) 

J-D 


were  evaluated  to  give  the  heat  content  and  potential  energy  respectively.  In  Eq.  (8) 
Cp  is  the  specific  heat  for  sea  water,  a  the  coefficient  of  thermal  expansion  and  D  the 
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depth  of  the  thermistor  chain  at  mooring  B.  The  error  in  potential  energy  chetnge 
due  to  neglecting  salinity  effects  on  density  is  estimated  to  be  less  than  10%. 

Finally,  the  third  set  of  parameters  (Fig.  32a, d)  comprise  plots  of  barometric  pres¬ 
sure  BP,  wind  speed  WS,  air  temperature  AT,  sea-surface  temperature  ST  and 
estimates  of  sensible  and  latent  heat  fluxes  QS  and  QL  respectively.  The  derivation 
oi  QS  and  QL  involved  first  estimating  QS  from  a  bulk  aerodynamic  formula  and 
then,  estimating  QL  using  an  empirical  formula  for  the  Bowen  ratio,  B  =  QS/QL, 
dependent  on  air/sea  temperature  difference  and  wind  speed.  Further  details  of  the 
method  are  presented  in  Moen  (1984).  Both  the  second  and  third  set  of  parameters 
(Figs.  31  and  32)  derive  from  instruments  at  mooring  B  since  these  presented  the 
longest  time  series. 

We  can  now  make  a  number  of  comments  regarding  the  parameters  plotted.  The 
Richardson  number  Ri  shown  in  Figs.  27  to  30  measures  dynamic  stability  in  the 
presence  of  stratification  and  is  thus  an  indication  of  mixing.  As  defined  by  Eq.  (7), 
Ri  is  a  bulk  value  over  a  10  m  or,  in  some  cases,  a  20  m  layer,  so  that  according  to 
PuUard  et  al.  ( 1 973)  mixing  will  occur  if  its  value  drops  below  unity  -  the  mechanism 
being  one  of  high  levels  of  production  of  turbulent  mixing  energy,  attributable  not 
only  to  the  mean  shear  but  also  to  breaking  internal  waves. 

Scanning  the  plots  (Figs.  27  to  30)  we  see  that  with  few  exceptions  Ri  remains 
below  1  both  within  the  mixed  layer,  as  expected,  and  across  the  mixed  layer  base 
during  the  entire  period  of  measurements.  There  are  also  periods  of  potential  mixing 
showing  up  in  the  thermocline.  For  example,  for  the  layer  40-60  m,  Ri  is  less 
than  1  more  or  less  continuously  during  days  285-287  and  periodicaUy  (with  inertial 
frequency)  during  days  290-294  and  297-300.  Similarly,  for  the  60-80  m  layer  Ri 
falls  below  1  on  three  occasions  between  days  285  and  287  and  for  about  30%  of 
the  time  during  days  2970-301.  Reference  back  to  Fig.  25  shows  that,  at  least 
for  the  period  days  285-287,  the  40  m  and  60  m  current  vectors  are  180°  out  of 
phase.  Verticcd  variations  in  phase  of  the  inertial  osciUations,  therefore,  can  have  a 
significant  effect  on  mixing  within  the  thermocline. 

It  is  perhaps  a  point  of  curiosity  that  with  the  bulk  gradient  Richardson  number 
being  so  frequently  below  unity  across  the  mixed-layer/ thermocline  interface,  more 
deepening  of  the  mixed  layer  did  not  occur. 

The  expression  for  changes  in  potential  energy  given  by  Eq.  (8)  assiunes  no  horizontal 
or  vertical  advection  is  taking  place.  The  expression  is  equivalent,  to  the  lowest  order 
in  aAT,  to  the  more  accurate  and  physically  meaningful  method  of  calculating 
potential  energy  where  vertical  integration  begins  at  a  geopotential  surface  and 
continues  up  to  the  free  surface,  thereby  allowing  for  expansion  of  the  water  column. 
Since  a  downwelling,  such  as  occurs  during  half  of  an  inertial  cycle,  will  cause  the 
measured  portion  of  the  water  column  to  be  warmer,  Eq.  (8)  wiU  yield,  incorrectly, 
an  increase  in  potential  energy.  That  is,  an  expansion  of  the  water  column  (through 
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heating)  is  assumed  to  have  taken  place.  Inertial  cycles  in  the  potential  energy  plots 
of  Fig.  31,  therefore,  are  meaningless  as  would  be  any  other  advective  effects  on 
potential  energy  estimates  by  Eq.  (8). 

A  possible  advective  effect  might  explain  why  for  the  period  days  294-296  the  heat 
content  and  potential  energy  both  increase  (Fig.  31)  despite  the  mixed  layer  having 
deepened  from  36  to  44  m  and  the  air  temperature  having  dropped  from  22  to 
19  °C  (Fig.  32).  By  day  300,  however,  the  heat  content  has  dropped  to  its  day- 
294  value  but  the  potential  energy  remains  high,  reflecting  the  mixing  that  has 
occurred.  Another  example  of  mixing  increasing  the  potential  energy  can  be  seen 
for  the  period,  days  280-285  (Fig.  31).  Here  the  mixed  layer  deepens  from  32  to  38  m 
but  the  maximum  buoyancy  and  heat  --ontent  show  only  slight  changes,  whereas  the 
potential  energy  increases  quite  noticeably. 

Returning  to  the  heating  of  the  water  mentioned  above  for  days  294-296  (Fig.  31c), 
it  should  be  pointed  out  that  despite  the  strong  cooling  by  latent  and  sensible  heat 
lossec  (Fig.  32c)  there  might  have  been  a  net  radiative  heating  for  the  same  period. 
There  is  no  direct  evidence  for  this  since  cloud  cover  statistics  were  not  taken,  but  one 
can  imagine  the  colder  and  perhaps  drier  air  to  have  been  associated  with  cloud-free 
skies.  Radiative  heating,  therefore,  presents  cm  alternative  speculation  to  advective 
effects  for  the  observed  increase  in  heat  content  and  consequent  increase  in  potential 
energy  for  the  two-day  period. 
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7.  Summary  and  conclusions 

In  this  section,  points  of  summary,  conclusion  or  recommendation  are  presented  in 
an  order  approximately  following  their  discussion  in  the  report. 


7.1.  SPATIAL  DESCRIPTION 

(1)  Two  hydrographic  surveys  were  conducted:  one  during  mooring  deployment 

43  CTD  stations  and  180  XBT  casts)  and  one  by  the  Italian  Navy  during 
late  October  (~  16  CTD  stations  and  120  XBT  casts.  These  data  form  the 
basis  of  the  spatial  descriptions. 

(2)  In  describing  the  spatial  structure  a  distinction  was  made  between  an  upper 
or  ‘surface’  layer  of  about  150  m  depth  and  a  deeper  ‘intermediate’  layer  from 
about  150  to  1000  m.  This  division  was  somewhat  arbitrary  and  although 
corresponding  to  about  the  lower  portion  of  the  seasonal  thermocline  it  has 
no  special  dynamicaj  significance.  Rather,  it  had  to  do  with  display;  the  very 
strong  vertical  gradients  of  the  surface  layer  required  large  contour  intervals, 
and  these  would  have  been  masked  if  the  layers  had  not  been  viewed  sepa¬ 
rately;  the  relatively  weaker  structure  is  present  in  the  intermediate  layer. 

(3)  The  upper  layer  showed  an  essentially  uniform  spatial  structure  in  the  study 
sone  with  no  large-scale  variations  in  the  south-north  direction.  However, 
shallowing  of  the  potentially  mixed  layer  was  seen  in  both  the  tempera¬ 
ture  emd  salinity  west-east  sections  on  approaching  the  Sardinian  shelf  break. 
Since  the  surface  temperature  was  also  warmer  in  this  region,  the  shallowing 
of  the  layer  (probably  caused  by  upwelling)  must  have  been  happening  slowly 
compared  with  the  heating  time-scale  through  insolation. 

(4)  The  general  pattern  of  geostrophic  flow  through  the  study  area  relative  to 
motion  at  1000  m  appeared  consistent  with  the  maps  presented  in  Ovchin¬ 
nikov  (1966).  In  the  upper  200  m  the  eastward  flowing  North  African  current 
branched  northward  into  a  portion  of  the  study  area  and  then  flowed  south¬ 
eastward  again  causing  an  anticyclonic  motion  which  appeared  to  have  almost 
pinched  off,  into  an  eddy.  At  the  20  m  level  a  partially-formed  cyclonic  gyre 
was  identified  in  the  northeast  portion  of  the  study  area.  This  was  associated 
with  the  cold,  possibly  upwelled,  water  showing  up  in  the  temperature  data 
referred  to  in  the  preceding  paragraph  (item  3).  The  geostrophic  flow  near 
the  mooring  sites  was  southward  at  about  7  cm/s. 

(5)  The  deeper  layer  (150  to  1000  m)  was  characterized  by  the  so-called  ‘Levan¬ 
tine  Water’  lying  approximately  between  300  and  650  m.  This  water,  show¬ 
ing  general  agreement  with  transects  presented  in  the  Mediterranean  Atlas  of 
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Miller,  Tchernia  and  Chctrnock  (1970),  indicated  a  weak  north-south  gradient 
in  temperature  and  salinity  -  hoth  decreasing  to  the  south. 

The  most  striking  features  of  the  deeper  layer  were  the  warm  and  cold  tem¬ 
perature  ‘blobs’,  density  -  compensated  by  their  salinity,  which  were  observed 
just  above  the  Levantine  Water.  Within  the  Levantine  Water  of  the  west¬ 
ern  portion  of  the  study  area  several  of  the  CTD  profiles  indicated  a  fairly 
pronounced  step-like  structure. 


7.2.  TEMPORAL  DESCRIPTION  (THERMISTOR-CHAIN  DATA) 

(1)  Time  series  for  the  approximately  40-day  period  were  collected  from  3  moored 
thermistor  chains,  each  chain  consisting  of  22  thermistors  spanning  mixed 
layer  and  thermocline  (to  ~  100  m  depth)  and  sepsu-ated  one  chain  from  the 
other  by  about  10  km,  in  a  triangular  arrangement. 

(2)  Periodic  vertical  displacement  of  isotherms  indicated  inertial  or  near-inertial 
frequency  motion  to  have  been  present  throughout  the  measurement  period. 
Indeed,  these  motions  dominated  the  records  -  particularly,  as  expected,  after 
the  onset  of  the  storm  events. 

(3)  The  amplitude  of  the  oscillations  was  frequently  seen  to  be  greater  in  the  more 
weakly  stratified  water  below  the  thermocline  than  in  the  thermocline  itself. 
The  total  energy  of  the  vertical  oscillations,  however,  may  not  have  been 
increasing  with  depth  in  the  same  way,  since  the  buoyancy  restoring  forces 
due  to  stratification  were  a  higher  in  the  upper  than  in  the  lower  thermoclin° 

(4)  A  significant  delay  between  onset  of  vertical  oscillations  and  sudden  increase 
in  wind  speed  (in  one  ceae  ~  24  h)  was  noted  on  several  occasions.  This 
observation  was  brought  to  bear  on  the  possible  generation  mechanism  for 
the  vertical  oscillations.  Four  possible  mechanisms  were  discussed: 

(i)  Advection  of  sloping  isopycnal  surfaces  by  horizontal  inertial  oscilla¬ 
tions. 

(ii)  A  horizontal  Ekman  mass-flux  divergence  due  to  a  spatially  varying 
wind  stress.  The  inertial  frequency  component  of  this  diverging  mass 
flux  then  forces  complementary  vertical  oscillations  (KroU,  1975). 

(iii)  A  parametric  instability  (Kroll,  1982  and  1984)  of  a  slab-like  mixed- 
layer  flow  (inertially  oscillating  or  not)  when  subjected  suddenly  to  a 
spatially  uniform  wind  stress. 

(iv)  Vertical  inertial  oscillations  generated  elsewhere,  particularly  by  mech¬ 
anism  (ii)  above,  may  also  have  propagated  into  the  study  area. 

Of  these  mechanisms,  (iii)  and  (iv)  were  the  only  ones  to  offer  explanation 
of  the  apparently  delayed  response  of  the  vertical  inertial  oscillations  to  the 
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sudden  strong  changes  in  the  wind  stress,  such  as  occurred  several  times 
during  MILEX-82. 

(5)  A  fifth  mechanism  capable  of  generating  vertical  inertial  oscillations  locally 
might  arise  from  step  changes  in  a  spatially  uniform  wind  stress  over  water 
with  a  high  background  vorticity.  This  might  occur,  for  example,  when  a 
storm  passes  over  a  geostrophic  eddy.  Such  a  mechanism  was  not  discussed 
in  the  text  of  the  report  and  was  unlikely  to  have  been  relevant  to  the  study 
area;  no  eddies  were  observed  either  by  satellite  or  from  the  hydrographic 
data  during  MILEX-82.  Details  of  the  potential  mechanism  do  not  appear  in 
the  literature  but  they  can  be  deduced  (not  to  be  presented  here)  from  the 
Ekman-layer/geostrophic-layer  modelling  presented  in  Stern  (1975). 

(6)  Vertical  oscillations  super-inertial  in  frequency  (u)  =  2/)  were  aJso  observed 
in  the  records.  As  has  been  pointed  out,  for  example  by  Krauss  (1981), 
Rubenstein  (1983)  and  KroU  (1984),  these  higher-frequency  components  are 
associated  with  high  downward  vertical  group  velocities  and  are  thus  efficient 
in  extracting  energy  from  the  mixed-layer  flow.  This  would  explain  why  the 
mixed  layer  did  not  deepen  more,  considering  the  severity  of  the  storm  events 
during  MILEX-82.  That  is  to  say  mixing  energy  may  have  been  propagating 
out  of  the  mixed  layer,  leaving  reduced  amounts  available  for  mixing. 

According  to  KroU  (1984)  super-inertial  waves  are  a  prediction  of  the  para¬ 
metric  instability  mechanism  (item  4  (iii)  above)  for  generation  of  vertical 
inertial  osciUations.  Their  presence,  therefore,  is  suportive  (although  not 
exclusively  so)  of  the  instability  mechanism. 

(7)  The  vertical  near-inertial  oscillations  did  not  display  a  constant  phase  with 
depth  and  the  phase  shifts  were  not  constant  with  time.  Although  a  180° 
phase  shift  is  predicted  smd  observed  (Krauss,  1979)  in  a  two-layer  shtdlow 
water  model  for  the  horizontal  osciUations,  it  is  not  clear  why  phase  changes 
with  depth  should  be  taking  place  in  the  case  of  deep  stratified  water.  Further 
and  more  detailed  anedyses  of  these  data  using  techniques  such  as  complex 
demodulations  are  recommended. 
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7.3.  TEWPO^AI.  LiESCRIPTION  (CURRENT-METER  DATA) 

(1)  Vertical  arrays  of  current  meters,  distributed  between  the  three  mooring  lo¬ 
cations,  together  provided  flow  information  for  the  upper  (150  m)  layer. 

(2)  The  mean  flows  are  relatively  weak  (in  the  1  to  10  cm/s  range)  at  all  depths, 
but  the  inertial  oscillations  were  of  quite  high  amplitude  (up  to  50  cm/s], 
particularly  in  the  surface  mixed-layer  after  sudden  large  increases  in  wlpfl 
stress. 

(3)  The  mixed-layer  flow  was  ‘slab-like’,  but  with  the  inertial  oscillations  at  the 
layer  base  (~  30  m)  mostly  out  of  phase  with  the  layer  itself  by  about  90’’  ± 
40°. 

(4)  Unlike  the  vertical,  the  horizontal  oscillations  responded  rapidly  to  the  sudden 
changes  in  wind  stress  of  day  279.  This  was  seen  in  the  mixed  layer  as  not 
an  amplitude  response  but  a  sudden  phase  change  of  the  existing  inertial 
oscillations. 

(5)  The  amplitude  response  was  found  to  be  slower,  with  a  delay  of  about  14  h. 
Two  reasons  for  this  are  possible:  firstly,  before  they  could  grow  in  zunplitude 
the  new  oscillations  had  to  overcome  or  cancel  the  existing  out-of-phase  os¬ 
cillations;  and  secondly,  there  was  a  marked  tendency  for  the  wind  vector  to 
rotate  in  a  counter-clockwise  sense,  thus  opposing  the  very  inertial  motions 
it  was  generating. 

(6)  An  interesting  feature  of  the  inertial  currents  was  in  the  variability  of  their 
phase  with  both  depth  and  time.  Sustained  periods  of  large  phase  diflTerences 
and  even  periods  of  antiphase  occurred  within  the  thermocline,  particularly 
in  the  40  to  60  m  depth  range.  The  resulting,  relatively  large  shear  caused 
lower  values  of  the  Richardson  number  in  this  region,  creating  the  conditions 
for  a  raised  level  of  turbulent  diffusion. 

(7)  The  variability  in  phase  was  itself  seen  to  be  largely  a  manifestation  of  changes 
in  the  near-inertial  frequency  with  depth  and  time.  In  the  mean,  there  was 
a  tendency  for  frequency  to  increase  with  depth  -  the  mixed-layer  value  of 
about  1.01/  being  less  than  the  pycnocline  value  of  about  1.05/.  Over  a  25- 
day  period  the  frequency  at  40  m  was  to  decrease  from  1.15/  to  1.0/  while 
at  60  m  it  increased  during  the  same  period  from  1.05/  to  1.15/. 

(8)  The  observed  value  of  the  near-inertial  frequency  in  the  mixed  layer  can  only 
differ  significantly  from  the  local  value  of  /  if  there  is  a  high  background 
geostrophic  vorticity  u/  present  in  the  water  column.  The  relatively  small  de¬ 
viation  of  the  near-inertial  from  the  inertied  frequency  observed  in  the  mixed 
layer  of  the  MILEX-82  study  suea  was  therefore  testimony  to  the  very  low 
background  or  geosptrophic  vorticity  present  there. 

(9)  Reasons  for  the  variability  in  near-inertial  frequencies  were  presented  as: 

(i)  the  tendency  for  higher  frequency  inertial  energy  to  have  faster  down- 
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ward  group  velocity,  causing  it  to  disappear  from  the  mixed  layer  more 
quickly;  and 

(ii)  deeper  inertial-energy  may  have  been  generated  well  north  of  the  study 
area  (as  mentioned  for  the  vertical  oscillations  -  Sects.  4.3  and  7.2) 
where  the  locsd  value  of  /  was  much  higher. 

(10)  Super-inertial  frequency  motions  were  seen  as  secondary  peaks  in  the  velocity 
spectra.  These  were  most  likely  a  result  of  strong  non-linearity  due  to  inter¬ 
action  between  the  vertical  and  horizontal  inertial  oscillations  (i.e.  terms  in 
the  momentum  balance  such  as  w  du/dz).  Estimates  from  the  data  showed 
this  term  in  the  mixed-layer/thermocline  interface  region  to  be  of  a  similar 
order  effect,  to  the  Coriolis  term  (fv)  itself. 


7.4.  OTHER  PARAMETERS 

( 1 )  Time-series  plots  were  presented  for  the  following  measured  or  derived  param¬ 
eters:  vertical  current  shear,  Brimt-Vaisala  frequency,  bulk  gradient  Richard¬ 
son  number,  mixed-layer  depth,  maximum  Brunt-Vaisala  frequency,  change  in 
water  column  heat  content,  change  in  water-column  potential  energy,  baro¬ 
metric  pressure,  wind  speed,  air  temperature,  sea-surface  temperature  and 
finally  estimated  sensible  and  latent  heat  fluxes. 

(2)  The  estimated  gradient  Richardson  numbers  for  2  levels  within  the  thermo- 
cline  were  found  to  be  below  unity  during  a  number  of  prolonged  periods. 
These  periods  coresponded  to  times  of  large  vertical  changes  in  the  phase 
of  the  inertial  oscillations,  confirming  the  importance  of  inertial-wave  phase 
variations  with  depth  to  mixing  within  the  thermocline. 
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Fig.  la.  Location  of  MILEX'82  CTD  stations  (large  numbers)  and  XBT  stations  (small 
numbers). 
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Fig.  I  b.  Location  of  MILBX-82  moorings. 
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Fig.  2.  Schematic  i«pre’<cutstion  of  the  six  moorings. 


SaCLANTCEN  SM-213 


SaCLANTCBN  SM-213 


5*  6®  7®  8®  9® 


Fig.  4a.  CTD/XBT  survey  track  of  the  Italian  Navy  hydrographic  vessel  Ammiraglio  Mag- 


Fig.  4b.  Ship’s  track  {Maria  Paolina  G.)  for  search  and  recovery  of  drifting  moorings. 


k 


F.  RESPONSE  FNC. 


Fi^.  5.  Plots  of  impulse,  response  and  frequency  response  for  the  wrapped  thermistor. 
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Fig.  6.  Vertical  temperature  sections  for  the  upper  150  m:  (a)  south-north  and  (b)  west 
east. 
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Fig.  14.  Example  of  stei>-ljke  structure  in  the  Levantine  layer  observed  for  CTD  stations  8 
and  30.  The  enlargement  of  steps  shown  in  potential  density  are  for  station  8. 
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Fig.  J5a.  Temperature  contours  from  thermistor  chain  A1  for  days  270-280.  Contour 
intervab  are  every  0.5  *C  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  li.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne. 
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Fig.  15b.  Temperature  contours  from  thermistor  chain  A  I  for  days  28H-290.  Contour 
interval  are  every  0.5  *C  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne. 
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Fig.  J5c.  Temperiture  contours  from  thermistor  chain  A1  for  days  290-300.  Contour 
intervals  are  very  0.5  °C  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocline. 
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Fig.  16a.  Temperature  contours  from  thermistor  chain  Bl  for  days  270-280.  Contour 
intervals  are  every  0.5  and  positions  with  temperatures  of  uppermost^and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne. 
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Fig.  16b.  Temperature  contours  from  thermistor  chain  Bl  for  days  280-290.  Contour 
intervals  are  every  0.5  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne- 
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Fig.  16c.  Temperature  contours  from  thermistor  chain  Bl  for  days  290-300.  Contour 
intervals  are  every  0.5  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne. 
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Fig.  J6d.  Temperature  contours  from  thermistor  chain  Bl  for  days  300-310.  Contour 
intervals  are  every  0.5  and  positions  with  temperatures  of  uppermost-and  lowermost 
thermistors  are  indicated  by  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocUne. 
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Fig.  17a.  Temperature  contours  from  thermistor  chain  Cl  for  days  270-280.  Contour 
intervals  are  every  0.5  ^C  and  positions  with  temperatures  for  uppermost- and  lowermost 
thermistors  are  indicated  with  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocline. 
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Fig.  17b.  Temperature  contours  from  thermistor  chain  Cl  for  days  280-290.  Contour 
intervals  ate  very  0.5  ”C  and  positions  with  temperatures  for  uppermost-and  lowermost 
thermistors  are  indicated  with  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocline. 
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Fig.  17c.  Temperature  contours  from  thermistor  chain  Cl  for  days  290-300.  Contour 
intervals  are  every  0.5  *C  and  positions  with  temperature  for  uppermost-and  lowermost 
thermistors  are  indicated  with  a  star  symbol  every  24  h.  Wind  speed  is  plotted  to  indicate 
wind  frequency  on  the  thermocline. 
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Fig.  J7d.  Temperature  contours  from  thermistor  chain  Cl  for  days  300-310.  Contour 
intervals  are  every  0.5  *C  and  positions  with  temperatures  for  uppermost-and  lowermost 
thermistors  are  indicated  with  a  star  symbol  every  24  h  Wind  speed  is  plotted  to  indicate 
wind  forcing  on  the  thermocline. 
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Fig.  I8a.  Temperature  variance  spectra  for  the  upper  12  thermistors  at  mooring 
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Fig.  1 9b.  Temp«iaiui«  variance  spectra  for  the  lower  10  thermistors  and  the  chain  average 
at  mooring  Bl. 
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Fig.  20a.  Temperature  variance  spectra  for  the  upper  12  thermistors  at  mooring  Cl. 
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Fig.  20b.  Tempetaiute  variance  for  the  lower  10  thermistors  and  chain  average  of  moor- 


(M)**2/CYC/HR 


SaCLANTCEN  SM-213 


WAVE  PER'lOD  (HRS) 

lo’  lo’  lo'  10°  lO'' 


FREQUENCY  (CYC/HR)  FREQUENCY  (CYC/HR) 


Fig.  21.  Spectta  for  isotherm  displacement  at  mooring  Cl. 
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Fig.  23a.  Composite  stick  diagram  representing  period  JD  270-280  for  all  eight  current- 
meters  with  depth  and  current-meter  type  indicated.  The  scale  is  indicated  at  the  bottom 
of  the  plot.  Wind  sticks  for  the  same  period  are  plotted  across  the  top  to  a  different  scale. 
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Fig.  23b.  Composite  stick  diagram  representing  period  JD  2b0-290  for  current  and  wind 
vectors. 


Fig.  23c.  Composite  stick  diagram  representing  period  ID  290-300  for  current  and  wind 
vectors. 
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Rotary  spectra  for  wind  energy  at  mooring  B 
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Fig.  25.  Unit  vectors  at  25-  intervals  comparing  current  direction  for  each  current  meter 
relative  to  the  direction  at  40  m.  Directions  taken  from  Fig.  23;  current-meter  depth  is 
written  on  each  stick. 


-  @33  - 


Fig.  27a.  Time  series  plots  representing  period  JD  270-280  of  bulk  estimates  of  Brunt- 
Vaisala  frequency  N,  shear  5,  and  gradient  Richardson  numb\.r  Ri  for  depths  10-20  in, 
20-30  m  and  30-40  m.  The  mooring  or  moorings  involved  for  each  series  is  indicated, 
preceding  the  specified  depth  range. 
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Fig.  28b.  Brunt-Vaisala  frequency  N,  shear  S  and  gradient  Richardson  number  Ri  for 
period  JD  280-290  and  depths  40-60  m  and  60-80  m. 
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Fig.  29.  Time  series  plots  representing  period  JD  290-300  of  bulk  estimates  of  Brunt- 
Vaisala  frequency  N,  shear  5,  and  gradient  Richardson  number  Ri  for  depth  40-60  in 
and  60-80  m.  The  mooring  or  mootings  involved  for  each  series  is  indicated,  preceding 
the  specified  depth  range. 


Fig.  30.  Time  series  representing  period  JD  300-310  of  bulk  estimates  of  Bruiit-Vaisala 
frequency  N,  shear  5,  and  gradient  Richardson  number  Ri  for  depths  40-60  m  and  SO¬ 
SO  m.  The  mooring  or  moorings  involved  for  each  series  is  indicated,  preceding  the  specified 
depth  range. 
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Fig.  31a.  Time  series  plots  for  period  JD  270-280  of  mixed-layer  depth  MLD,  maxi¬ 
mum  Brunt-Vaisala  frequency  ATmax,  change  in  heal  content  HTt-HTo,  and  change  in 
potential  energy  PEi~PEo.  All  estimates  are  for  mooring  Bl- 


Fig.  31b.  Time  series  plots  for  period  JD  280-290  of  mixed-layer  depth  MLD^  maxi¬ 
mum  Brunt-Vaisala  frequency  TVmax,  change  in  heat  content  HTX-HTo,  and  change  in 
potential  energy  PEt-PEo.  All  estimates  are  for  mooring  Bl. 
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Fig.  31c.  Time  series  plots  for  period  JD  290-300  of  mixed-layer  depth  M LD,  maximum 
Brunt-Vais^a  frequency  yVmax,  change  in  heat  content  HTx~HTo,  and  change  in  poten¬ 
tial  energy  PEi-PEo.  AU  estimates  are  for  mooring  Bl. 


Fig  3Id.  Time  series  plots  for  period  JD  300-310  of  mixed-layer  depth  MLD,  maxi¬ 
mum  Brunt-Vaisala  frequency  Nma.x,  change  in  heat  content  HTt-HTo.  and  change  in 
potential  energy  PEt-PEo.  All  estimates  ate  for  mooring  Bl. 
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Fig  Ti.iie  series  plots  for  period  JD  270  280  of  barotnt'.ric  nressure  BP.  wind  speec 

W'5,  air  lem  peralure  AT,  sea  surface  temperature  .S  I .  and  estimates  of  sensible  and  latent 
heat  fluxes  QS  and  QL  respectively-  All  estiniafe«;  are  f-T  nn'oriiifi  Rl 
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Fig.  32b.  Time  series  plots  for  period  JU  280  290  of  barometric  pressure  HP.  wind  speed 
IV  S,  air  temperature  AT,  sea  surface  temperalure  ST.  and  estimates  of  sensible  and  l-.ienl 
lieal  fluxes  QS  and  Qi  respectively  All  esliinales  are  for  nuiotini:  Hi 
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